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I consider an extented version of Bekenstein's Tensor- Vector-Scalar theory where the action of the 
vector field is of a general Einstein-Ether form. This work presents the cosmological equations of 
this theory, both at the background and perturbed level, for scalar, vector and tensor perturbation 
modes. By solving the background equations in the radiation era analytically, to an excellent 
approximation, I construct the primordial adiabatic perturbation for a general family of scalar field 
kinetic functions. 
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' In the last two decades, cosmology has undergone a "precision" revolution, with a large influx of data such as 
^ ^ , observations of the Cosmic Microwave Background large scale structure [3], and supernovae observations 0]. 

There is now a consensus cosmological standard model, named ACDM, based on General Relativity as the theory of 
\ gravity, which requires only about 4% of the energy budget of the universe to be in known baryonic form, while the 
' ■ rest is divided into two apparently distinct, dark components: cold dark matter and dark energy. It is unfortunate 
I— i' however, that apart from their phenomenology as fluids, we know nothing of their actual nature at the present time. 

, Cold dark matter, typically composed of very massive slowly moving and weakly interacting particles, is required on 
Oh cosmological scales mainly to source large scale structure. To dramatise its importance, assuming general relativity, 
Q , if dark matter was absent, structure as we know it would not have even formed yet. A plethora of such particles 
5-H • generally arises in particle physics models beyond the standard model quite naturally with the right cross-sections to 
' create the right abundance (see [3,01 for a reviews). Yet, while its phenomenology as a dust fluid has been shown to 
agree with observations to a very good degree, the actual nature of cold dark matter is left to speculation as no cold 
dark matter particle has been observed so far. Moreover there are still some mishaps within the ACDM paradigm, 
y—{ . for example the problem of voids [1] and the recent observations of the Abel 520 cluster Q . 

^ ' Given that the law of gravity plays a key role, to all observations from which dark matter and dark energy are 
inferred, it is conceivable that General Relativity breaks down at small enough gradients and curvatures, and an 
alternative theory of gravity might also provide an explanation to the dark sector. One such theory proposed some 
time ago by Bekenstein [s], building on key work by Sanders P. This theory was dubbed Tensor- Vector-Scalar 
(TeVeS) because it relies on a bimetric transformation involving a scalar and a vector field. It was designed to reduce 
to the Aquadratic Lagrangian non-relativistic theory of Bekenstein and Milgrom [lOll. Thus, it provides essentially 
the same phenomenology as Milgrom's Modified Newtonian Dynamics (MOND) [llj for galactic rotation curves, for 
which MOND has had tremendous success p^ . 

TeVeS theory has since been shown |l3ll to be able to source structure in a similar way as dark matter. The vector 
field in the theory plays a key role [lj| , as for a wide range of parameters it has a power-law growing mode which 
sources potential wells. In contrast with dark matter, the vector field has shear which creates a mismatch between 
' the two scalar gravitational potentials. This has been identified as prospective discriminator between dark matter 
I and theories like TeVeS [Tsl [l^ . [TtI. fisj. Other non-cosmological tests of TeVeS have also been studied, for example 
gravitational lensing 1^ 2^. Probing the difference between the arrival times of neutrinos and gravitational waves 
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from distant supernovae is also a possibility [2l|. 

Having shown that one can cast TeVeS in a single metric form, with the scalar field absorbed into the vector 
field [22], Zlosnik, Fcrreira and Starkman, have explored a sister theory, based solely on a unit-timelike vector field 
with a non-canonical kinetic term (23| which is a non-canonical Einstein-ether theory ^A] . This theory has also been 
shown to source structure in a similar way as TeVeS [25| (cosmology within the context of canonical Einstein-ether 
theory has been extensively studied in [26l. [27l [28j ) while its predictions for corrections to Newtonian gravity in 
the solar system have also been studied [2^. Further exploration of this theory into different directions has been 
considered by Halle and Zhao [13, HH . 

In this work I initiate a study of a version of TeVeS which involves a generalization of the action of the vector field 
into the Einstein-ether form. This is motivated in part from the instability present in spherically symmetric solutions 
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to the original TeVeS theory [s?] , that possibly stems from the fact that the kinetic term for the unit-timelike vector 
field was of Maxwellian form which can violate the dominant energy condition [ssj . Seifert has shown that more general 
Einstein-ether actions can be stable depending on the parameters. It is therefore of importance to check whether 
forms of TeVeS exist which have stable spherically symmetric solution and still can form large scale structure. Possible 
directions for deriving TeVeS-type theories from more fundamental theories are discussed in[34. ,35, ,36]. 

In section-im I lay down the action and derive the field equations. The cosmological equations for a (possibly 
curved) Friedman-Lemaitre-Robertson- Walker (FLRW) metric are studied in section- HTTl where it is shown that they 
are identical with the TeVeS equations upto a rescaling of Hubble's constant. The cosmological solutions are therefore 
the same as those in TeVeS. I conclude that section by deriving an approximate (to an excellent degree) solution in 
the radiation era for a general family of scalar field functions. 

To study large scale structure we need the cosmological perturbation equations about an FLRW universe. The 
gauge form- invariant perturbed equations are derived using the techniques in [37| and are given in section- HVl for all 
types of perturbations, namely scalar vector and tensor perturbations. Using the approximate background solution 
in the radiation era from section- Hlll I construct the primordial adiabatic perturbation. The construction of the most 
general, regular primordial perturbation is quite involved and is given elsewhere [38i |. Throughout the paper I use the 
conventions of Wald ^ . 



II. FUNDAMENTALS : ACTION AND FIELD EQUATIONS 



A. Preliminaries 



TeVeS theory and the generalization herein, is a bimetric theory where gravity is mediated by a tensor field gab 
with associated metric-compatible connection Vq and well defined inverse g"''' such that g°''^gcb = f^^f,, a timelike 
(dual) vector field Aa such that g°'^AaAh = — 1, and a scalar field (p. Matter is required to obey the weak equivalence 
principle, which means that there is a metric gab with associated metric-compatible connection Va, universal to all 
matter fields, such that test particles follow its geodesies. The tensor field gab will be called the Einstein- Hilbert frame 
metric (see below) while gab the matter frame metric. 

The relation between the four above tensor fields (when the field equations are satisfied) is 

9ab = e^^'^gab - 2 sinh(20)A,Af, (1) 

with inverse 

gab ^ g20~a6 _^ 3 sinh(2(/)) A"^ A'' (2) 

where A" = g°'''Ab. 



B. The action principle 



The theory is based on an action S, which splits asS = Sg + SA + S^ + Sm, where Sg,SA,S^ and Sm are the actions 
for gab, vector field Aa, scalar field <j) and matter respectively. 

The action for gab, Aa and is most easily written in the Einstein-Hilbert frame, and is such that Sg is of Einstcin- 
Hilbert form 



(3) 



where g and R are the determinant and scalar curvature of g^jy respectively and G is the bare gravitational constant. 
Due to the complicated nature of the equations, the numerical value of G will not be the measured value of Newton's 
constant as measured on Earth. The precice relation between them depends on the spherically symmetric solution 
which apart from depending on the arbitrary function V (see below), is not expected to be unique, just like the case 
of standard TeVeS ji,!!^. 

The action for the vector field Aa is given by 
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i^°''^'*V,AbV,Ad - 2X{AaA^ + I] 
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where 

j^abcd ^ {g'^'g^'^ - g'"^~g'"') + K+ (g'^'.g'"^ + g'^'^g^") + K^g'^^'g'"^ + Ka^'^A^A" (5) 

A is a Lagrange multiplier ensuring the timelike constraint on Aa and Kb, Kj^, Kq and Ka are dimensionless constants. 
The action for the scalar field is given by 

S^ = - ^ J d\^g [^i {r" " A'^A') ^acp^bcp + Vifi)] (6) 

where /i is a non-dynamical dimensionless scalar field, and V{p) is an arbitrary function which must be such that 
^ ^ /i^ as /z ^ in order to have exact MOND limit, while it must diverge as /x ^ /xq where fiQ is a constant, in 
order to have exact Newtonian limit [4lj . One example is the form considered in [4l| which is 

dv _ ^il A' 



2 .(A -Ma)" (7) 



dfj, 16n£g /t — 1 

where £b is a scale, /ia is a constant, n is an integer power and jl = -j^. This general class of functions will also be 
used in this work. 

The matter is coupled only to the matter frame metric gat and thus its action is of the form Sm[g,X^] = 
/ d'^x \J—g L\g,x^\ for some generic collection of matter fields . 

One can further generalize the vector field action by making the constants Kb, K^, Kq and Ka functions of the 
scalar field (j>, as well as generalizing the scalar field action by inserting (/)-dependent functions as coefficients of the 
terms g"''' and A°'A'' in the kinetic term and a potential for . I leave this for future investigations (if warranted). 



C. The field equations 

Variation wrt the Lagrange multiplier A gives back the timelike constraint on the vector field, g'^^'AaAi, = —1. 
The matter stress-energy tensor Tab is defined by varying of the matter action w.r.t the matter frame metric as 
SSrn = -Ud''x^TabSg-'. 

Define tensors S^"^^^^^ and J^^^^^ as 
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5\J\)A'^A^ + ~9f'5\^A,)A^ + -gf'S^^A.^A^ 
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KAg'^ig'^'A" + g'^'A") 



(9) 



respectively. 

Then the field equations for ijab are given by 



Gab = SttG [Tab + 2(1 - e-4^)A=r,(,^)] - XAaAt + M [^a^yb^ - 2A=Vc0^(aVfc)0j + - i^iV' - V) gab 
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where Gab is the Einstein tensor of ijab- 

The field equations for the vector field Aa are 



^ jabce ja ceb 

2"^ d~ -J d 



VaAbVcAe - XAd - nA'Vb(bVd(l) + 87rG(l - e~^'f')A''Ti 



(11) 
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The field equation for tlie scalar field (j) is 



Va 



11 - A" A'') Vb(f> = STrGe-^"^ [g"'' + le''^'^ A'' A^\ T^b 
where the non-dynamical field jj, is found by inverting 

- A-A"-) Va<l>Vb<l> = -V' 



(12) 



(13) 



and therefore the arbitrary function V and its derivatives are nothing but functions of kinetic terms for (j), contracted 
with g"^ and A". 



III. FLRW COSMOLOGY 



A. Equations 

A most convenient coordinate system that is commonly used in cosmological perturbation theory is the conformal 
synchronous coordinate system with t denoting conformal time and a;" the spatial coordinates. This gives the matter 
frame metric with scale factor a as 



ds^ = [-dt^ + qijdx'dx^] 



(14) 



where Qab is the metric of a space of constant curvature ^ , with radius of curvature Vc and where K = Q for a flat , 
ii" = 1 for positively curved and iiT = —1 for negatively curved space. The scale factor of the Einstein- frame metric is 
h = ae'^. The vanishing of the Lie derivative with respect to all the Killing vectors of the background spacetime gives 
(j) — (j)(t) only, while the vector field is pure gauge for this background. 

The scalar field is governed by the TeVeS constraint which in this coordinate system reads, 



2 d/j. 

which must be inverted to get p,{a, (j), ^), and the second order equation 



(15) 
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3ji-(l) + ^TvGa^e-'^'t'ip + 3P) 



(16) 



where U = fi + 2^/^^. Both of the above are unchanged from TeVeS, because they are not affected by the vector 
field action. 



Defining the constant Kp = 1 + Kq + ^K^, the Friedmann equation gives 



while the Raychandhuri equation is 



SttGP 



K 



(17) 



(18) 



where p and P are the energy density and pressure of a matter fluid, and evolve as 



;^ + 3-(l + u;)p = 0. 



(19) 



with w = P I p. 



B. Solution in the radiation era 



For the form of the function above (O, it has been shown [H, [3)[1H that the scalar field tracks the dominant fluid 
such that the energy density of the fluid relative to the energy density in the scalar field is constant. This tracker 
behaviour is found in many scalar field dark energy (quintessence) models 42]. However it has also been shown 14 1 
for the special case of the Bekenstein toy model (n = 2 and /ia = 2 in the function above), that in the radiation 
era for realistic baryon and radiation densities [43| the radiation tracker is almost never reached until just before the 
transition to the matter era. Instead the solution is such that 4> evolves as a powerlaw of the scale factor. I show here 
that this is also true for the generalized function above. 

In the deep radiation era, for the function (l7|), /i is very large and we get that 

_ H 1 + n 

During this time (f> is very subdominant and we may assume that | ~ ^- The Friedmann equation then gives 

3Kf^ = SirGa'^e-^'t'^Pr (21) 
where is the initial condition of and pr is the radiation energy density. Define 

- ZKpHl ^^^^ 

where por is the proper radiation density today (as given by the radiation temperature), is the Hubble constant 
today, and the solution to the Friedman equation during radiation era is 

a = ^/n^rHQt (23) 

Although (j) is subdominant, we expect it to grow in order to approach the tracker solution. We therefore assume 
the ansatz 

= 0,- 0ia" (24) 

which gives 

| = -m<?!>ia"- (25) 
a 

where 0i and m are constants. I assume that initially (j) ~ Q. Now the variable q = —2/10 evolves as 



which gives the evolution of p. as 



The solution is 



so 01 > and we get 



dq a , , 

—^ + 3q = 6KF- (26) 
a ma a 



+ {m+l)p= —-Kpe-'^''"' (27) 



din a m0i 



/X = -^Kpa-"" (28) 
rrupi 



m{(pi - (p) 



irrespective of Ib and /iq. We now use the TeVeS constraint on for the large p. to get 
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Using the solution a{t) to get (j>{t) and fi{t) and then use the constraint we get 

4 



m 



and 
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For the TeVeS case, i.e. n 

1 1/5 
27 ' 



2, and 



l/(3+n) 

0, we get the Dodelson-Liguori [3| result m 



(31) 



(32) 



4/5 and 



IV. COSMOLOGICAL PERTURBATIONS 



The perturbed equations are written directly in Fourier space, and are thus dependent on the wave number k. 
Let us first define the following constants : 

Kt = Kb + K+ - Ka (33) 
- K+ + ]^Ko (34) 
Kf = l + K^ + Kd (35) 
Rk = 1-^-^ (36) 

As it turns out, for scalar and tensor modes Ka never appears by itself to linear order and is always absorbed in to the 
combination Kt defined above. Moreover the constant Hd functions as a damping constant in the scalar mode vector 
field equation when expressed in gauge invariant variables. When Kd = Q, the vector field equation is independent of 
k. Thus Kd must be positive and very close to zero, for if it were negative, the perturbations would have negative 
square speed of sound on fiat space and thus be greatly unstable on small scales, while if it were large and positive, it 
would damp the vector field on cosmological scales, which would render the theory irrelevant for structure formation. 
The constants Kt and Rk must be obey Q < Kt < 2 and < Rk < 1 respectively for the energy densities and the 
square speed of sound of all modes to be positive (the last conditions are sufficient to ensure that > 0). Thus 
scalar modes depend on Kt, Kp and Rk and tensor modes on Kp and Rk only. 

For vector modes the situation is slightly different and all constants Ki are needed. In that case one can use the 
parameters Kt, Kp and Rk defined above but one more is needed which can be either Ka or Kb which must obey 

^ 2KfRk ■ 



A. Scalar modes 



Scalar modes are defined as in [37|. The scalar field is perturbed as = + (p. The vector field has only one 
scalar mode, a, because the other one is fixed by the timelike constraint and is defined as Ai — ae^'^Wia. The 
matter frame metric has four scalar modes S, C and v such that 500 — — a^(l — 2S), goi — ~ct^ViC: 9ij = 
a^{l + ^x)Qij + ^^{(li^^/ ^ \(lij<l^^y^ iv- The Einstein frame metric has also four scalar modes S, x, C, and v 
defined in a similar way, related to the matter frame metric asS^S + ^s, x = x + C = C ~ (1 ^ e~'^'^)a and 
u = V. The fluid variables are the density contrast 5 = ^ and momentum divergence 6 such that the fluid velocity 
perturbation is deflned as Ui = oSI iO, where Ua is the unit-timelike wrt gab fluid velocity. 



1. Fluid equations 



The density contrast evolves as 



<5 = 3^ (li; - Cl) 5+{l + w) (-k^e -\x + k\ ) (37) 



where = ^ is the speed of sound. 
The momentum divergence evolves as 



' l + w l + w 3 
where S is the fluid's scalar anisotropic stress (see (37|). 

2. Scalar field equations 



The scalar field perturbation evolves as 

ip = ^7 - 0s 

where the auxiliary scalar field perturbation 7 (4^ evolves as 



(1 + 3C^) J +4 + 8nGa^e-^^-^{p + 3P) 
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{1 + '3C^)S ~ (1 + 3w)iE + 2ip) 



3. Vector equation 



The vector field equation is 



where the auxiliary gauge invariant vector mode E evolves as 



E+-E 
b 



K 
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a. 



+ ^ ( f - Kpe"^ ) (y + 2(C - a)) - e^'^na ( k 



b 9 



X + k^v + 6-S + 6 -- + 2— - 2-0 a 



-f (26^^"^ - f )/i0 - 0a) 



= 



The coupling to the matter velocity in the equation above has been eliminated with the use of (jH]) below. 

^. Einstein equations 



The two Einstein constraint equations are 

b 



-Kp 



b~ 



X + 2fc^(a-C)+6-S 



SirGa^e-^'I'p \S - 2ip] + e-^'''k^(KtE + 2-a) - 7 





and 



{i + ei^)+2lE 
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SnGa^e-^'^ {p + P)0 + 2fL(j)ip j-e-'^'^a 
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while the propagation equations are 



6i 



Kf< X + 2k'ia - + 6-E + 6 



2--- 



ie-4^(fc2 _ ^)(^ ^_ ^2^) ^ 3^^^ ^ 2ATTGa^e~^^piC^J - 2w^) - 2k^e-''^E ~ 2k''e-^^\a (45) 



for the couphng to the perturbed pressure and 



i> + 2(C - d) + 2 U + - (;> + 2C - 2a) 
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for the couphng to the shear 
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B. Vector modes 

Vector modes are defined as in [s?!. All vector modes have two polarizations and are purely spacial and divergenless. 
The vector field has a vector mode Pi and is defined as Ai = ae~'^Pi. The matter frame metric has two vector modes 
Ti, and fi such that goi = —a?ri and gij — 2o?'S/ i^.ifjy The Einstein frame metric has also two vector modes fi, and fi 

defined in a similar way, related to the matter frame metric as fi — Vi — {\ — e^*'^)/?^ and /; — fi. The fiuid variable 
is the vector mode Vi in the fiuid momentum such that the fluid velocity perturbation is defined as Ui — avi. 

1. Fluid equations 



The fiuid vector mode v evolves as 

, a w 

l-3u; - + - 

^ 'a l + w 

where cr*^"-' is the fluid's vector anisotropic stress (see (37|). 



k^ - a(^) 



(47) 



2. Vector field equations 



The vector mode equation is 



while the auxiliary vector mode e evolves as 



f3 = e+[c^~-](3 



b . &2 h 



(48) 



-K,{e+ le] - 1 [1 + i2KB - l)e'^*] (k^ + ^y+l{l- KfRkc^^) ( k^ - ^ ] { f + r ~ f3 



+ [Kpe^^ - 1) y2- - 4^ + 4-^j /3 + {2e^^ - l) /x^^^ = 

Notice that unlike scalar modes, the parameter Ka is no longer redundant. In this case we can parametrize the vector 
field with Kt, Kb, Kp and Rk- 



3. Einstein field equations 



We have the constraint equations 



rf. 



KfRk (f + r-p)+ e-^^p] = -l&T:Ga^e-^^{p + P) 
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and the propagation equation 



KfRk 



f + r-f3 + 2 




f + r-(3 



(v) 



Tensor modes 



Tensor modes are defined as in 37]. All tensor modes have two polarizations and are purely spacial and divergenless. 
The matter frame metric has a tensor mode iJij,such that gij — a^Hij. The Einstein frame metric has also a tensor 
mode Hij defined in a similar way, related to the matter frame metric as Hij = Hij . The equation of motion for the 
tensor mode H is 



KfRk 



H + 2 




2K 



(T) 



where tr'-^-' is the fluid's tensor anisotropic stress (see [st}). 



D. Adiabatic initial conditions for scalar modes 



Conformal synchronous gauge in radiation era 



We start by adopting the scalar mode perturbation equations, to the synchronous gauge (defined as = = 0, 
X — h and —k'^v — h + Grj), in the radiation era, using the background solution discussed above. Lets also define 



the following dimensionless variables 



kt, V = k9, u — ka, a = -nk T, and y 



— 7 
P" 



All equations are then 



written in dimensionless form, where derivatives wrt x are denoted by a prime. For simplicity let us also define 
0,. = m(j>iilQl/^ Hjp where Hk — Ho/k, such that (f)' — —(j>rx"^~^ (see section- ITlI B[) . 
We have the fluid equations for photons given by 



-h' 



(49) 



and 



(50) 



where the photon shear as well as higher moments of the Boltzmann hierarchy are vanishingly small due to the 
tight-coupling of photons to baryons and arc ignored. Likewise the fluid equations for neutrinos are 



I"' 



(51) 



vl = ^Si, - (Ji, (52) 

and 

< = ^«. + ^(/i' + 6r?') (53) 

where higher moments of the Boltzmann hierarchy are small because they are of higher powers in expansions about 
X and are ignored. 

The scalar field evolves according to 

V' ^^^y + ^rx"'-'^ (54) 
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and 



y = 



(1 + 2C^) SKf 



X ' X X 



-y 



DrX 



2 - e-^^' 



The vector field obeys 



u = E — (f u 

X 



and 



1 



E' + -E 

X 



Kt 

+2{KFe^^' ~ 1) 



e^'I'^nd [h' + 2u] + 2(1 - KFRKe^'^')ri' 
3Kf 



, 1 2 

If H — tp H — 
X x"^ 



X 



(2e'*"^' - 1) {ip + (l>rx'^-^u) = 



Finally we need the two Einstein constraint equations 
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h' + Qip' + 2mH LP 

X 



iKr. 



and 



x^ X 2Kf 



(55) 



(56) 



2. Adiabatic ansatz 



The adiabatic mode is such that rj —t \ for a; — )■ while all other perturbations vanish in this limit (regularity 
assumption). The adiabatic mode ansatz -q — l+r]2X^ , h — h2X^ solves the matter equations to give 5y = 5^ = —^h2X^, 
v-^ = —jgh2X^, v„ — —jgh2X^ — ^(^2 + 6772)2;^ and Ui, — j^{h2 + &ri2)x'^ . We seek solutions to the scalar and vector 
field variables which are regular as x — > 0. All the scalar and vector field terms in the Einstein constraint equations 
are then subdominant to lowest order in x and we can solve them to get ho = — and m — t^/Tr p i To n ^2 where 



Now consider the scalar field equation ([5^ where the first term clearly dominates at early times because of the 
regularity condition on p and y. Let y ~ yox^ for some power p > to leading order. We can then solve ([51)1 to get 

and using the above solution along with the one already found for h and S into (I55p we get 

{1 + p + 2C4,)yoxP = 10KF(f>rh2x"' -SKpe-'^^^pQX^+P + 18{3+p)KF(t>rVox"'^'' 



+3KF(t)r 2-e 



ux""-' (58) 



Since all terms above save the last one are regular as x — s- 0, then the w-term must also be regular which means that 
u — uqx^ with Z + TO > 1. Since < to < 1 and p > , we have that 2 + p > to and m + p > m which means that as 
X — > 0, the first term would always dominate over the second and third term and the above equation is reduced to 

{1+P + 2C^)yoxP = \QKF(t>rh2x"' + mF<j)r [2 - e-*^^] uox'+'""i (59) 

Now, the vector field equations become 

{I + l)uox^-^ = Eox'i - pqx^+p (60) 

and 

Kt{l + q)Eox'i + 2e^'l'^Kdh2X^ +A{1- KFRKe'^^')mx^ + '^{KFe^^' -1)uqx^-^ = (61) 
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where I have ignored the tp-terms because they are all ~ and so the h and 77-terms always dominate them as 

X ^ 0. I have also kept only the leading u-term. Using ([SO]) into (pT|) we get 



Eqx'^ + 2e'^^^Kdh2X^ KFRKe^^')mx^ = 



(62) 



where once again the (f-teim. has been ignored as it is of higher order than the h and 77 terms. Therefore consistency 
requires that q = 2 from which we get that the (^-term is of higher order than the ii^-term in (|5D1) which gives / = 3. 
Hence, the u-term in (|59p is of higher order and we get p = m. 

Reconstructing the full solution by using the above powers and matching coefficients then gives 



h 



7 



2Kf 



10-15Rk-4S, 
6(15i?K + 45,) 



—ht 
18 

15Rk 



45, 



18(15i?K + 45',; 
4 



ht 



E = 



3(15i?K + 4S',) 

20{KpRKe'^^^ - 1) + 2{15Rk + 4.S^){1 - e^'^^Kp) 



3(15i?K + iS^){3Kt + Kpe^'^' - 1) 



1 



1 + m + 2Crf, 



6(2 + m)ifi? 



It' 



Note that for standard TeVeS with 0^ = we get that E — and a = to this order. In this very special case 
E = 0(3) and a = 0(4) and depend on higher powers of h and rj. 



V. CONCLUSION 



I have formulated the cosmological equations both at the background and linear perturbation level for a version of 
TeVeS theory with a generalized vector field action. Using an analytical solution to the background equations for a 
general family of scalar field functions I constructed the primordial adiabatic perturbation. The most general type of 
regular primordial perturnation is studied elsewhere [ssj . These equations can be used to study large scale structure 
for these theories, to check whether there are stable versions of TeVeS which can agree with observations. 
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